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Rare Decays of D Mesons 
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We analyze the dominant decay mechanisms in the radiative decays D — > Vy, D — > P(V)f + f~, D — > 77, D — > £ + £~y, and we discuss 
possible chances to observe the effects of new physics in these decays. Then we comment on Cabibbo allowed D — > fory decays with 
nonresonant for and we probe the role of light vector mesons in these decays. 
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The K and B meson rare decays offer a chance to investi- 
gate effects of new physics in the down-like quark sector, 
while the rare D decays offer an opportunity to investigate 
the FCNC effects in the up-like sector. For the rare D de- 
cays only upper bounds have been established so far for a 
sizable number of the radiative or dilepton weak decay of D 
mesons [ Q] |2| |3] ID . New chances to observe these decay 
rates will appear soon at BELLE, BABAR and Tevatron. 
In analysis [ [5] |6| [7] ID ID of the dominant long distance 
contributions in D -> Vy, D -> P(W)£ + £ ~, D -> yy, and 
D — > t + t~y were determined. Possibilities to observe the 
effects of new physics were discussed. 

In the standard model (SM) the c — > uy transition is 
strongly GIM suppressed at one loop electroweak order 
giving a branching ratio of order 10 18 [ 10 1 . TheQCD cor- 
rected effective Lagrangian gives BR(c — > uy) 3 x 10~ 8 [ 
II II . A variety of models beyond the standard model were 
investigated and it was found that the gluino exchange di- 
agrams [ 12 1 within general minimal supersymmetric SM 
(MSSM) give the largest enhancement 
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The c — > xx£ + €~ amplitude is given by the y and Z penguin 
diagrams and W box diagram at one-loop electroweak or- 
der in the SM, and is dominated by the light quark contri- 
butions in the loop. We have used renormalization group 
improved effective weak Lagrangian within SM [ 9 1 and we 
found that the rates are BR(c -> ue + e~) S M = 2.4 x 10~ 10 
and BR(c -> uju + /i")sm = 0.5 x 10~ 10 . However, the 
present bounds on MSSM parameters lead to BR(c — > 
ue + e~) M ssM = 6.0 x 10~ 8 and BR(c -> u/i + /T)mssm = 
2.0 x 10~ 8 [Q. 

Long distance physics overshadows such effects since it 
usually dominates the decay amplitude [ 5 1. The long dis- 
tance contribution is induced by the effective nonleptonic 



|Ac| = 1 weak Lagrangian 

L - - Aij[aiurf{\ - y 5 )q i q j y fl (l - y 5 )c 

+ (hq/YnO- - 75)#;«y(l - ys)c], (2) 

with Aij = {GfI V2) VlqVuqi, accompanied by the emis- 
sion of the virtual photon. Here q u j denote the d or s quark 
fields. The effective Wilson coefficients are a\ = 1.2 and 
«2 = -0.5 [|5]. In our calculations of the long distance 
effects we use [01610 tne theoretical framework of heavy 
meson chiral Lagrangian M13I . The factorization approx- 
imation has been used for the calculation of weak transi- 
tion elements. In Table ^ we gi ve the branching ratios 
of D — > Vy decays [[5). The uncertainty is due to rela- 
tive unknown phases of various contributions. Although 



Table 1. The branching ratios for D — » V7 decays. 



D -> Vy 



BR 





K*°y 


[6- 


36] x 10- 5 




^p + y 


[20- 


- 80] x 10~ 5 


D° 


-+p°y 


[0.1 


- 1] x 10~ 5 


D° 


— > cjy 


[0.1 


- 0.9] x 10~ 5 


D° 


-> <5y 


[0.4 


- 1.9] x 10~ 5 


D + 


-» p + j 


[0.4 


-6.3] x 10- 5 


D s + 


— > K* + y 


[1.2 


-5.1] x 10~ 5 


D + 


K* + y 


[0.3 


-4.4] x lO" 6 


D° 


^K*°y 


[0.3 


-2.0] x 10~ 6 
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the branching ratios are dominated by the long distance 
contributions, the size of the short distance contribution 
can be extracted from the difference of the decay widths 
r(D° -> p°y) and F(D° -> ojy) [ 14]. Namely, the long 
distance mechanism cu — > ddy screens the cu — > uuy tran- 
sition in D° — > p°y and D° — > a>y, the p° and <±> mesons 
being mixture of uu and dd. Fortunately, the LD contribu- 
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tions are mostly canceled in the ratio 

BR(D° -> p°y) - BR(D° 



R 



ury) 



oc Re 



BR(D° -> uyy) 
A(D° -> uuy) 
A(D° ddy)' 



(3) 



which is proportional to the SD amplitude A(D° — > utiy) 
driven by c — > uy. This ratio /?sm - 6 + 15% in the Ref. [ 
fill, and can be enhanced up to 0(1) in the MSSM. In addi- 
tion to the c — > uy searches in the charm meson decays, we 
have suggested to search for this transition in B c — > B*y 
decay N15I . where the long distance contribution is much 
smaller. 

The amplitudes for exclusive decays D — > W€ + (~ and 
D — > Pl+f are dominated by the long distance contribu- 
tions. The branching ratios for these decays, as obtained in 
[|6||7 1, are given in Tables[5]and[3] In these Tables the first 
column represents the SD contributions, while in the sec- 
ond column the rates coming from the LD contributions are 
given. The rates for D — > Ve + e~ are comparable to those 
in Tableland can be found in [|6|. 



Table 2. The D — » V// + // branching ratios. 



D -> V/i V 



BR 



SD 



BR 



LD 



D° 


-> K* ( V/ V 







[1.6 


- 1.9] x 10" 6 


D+ 


-* p W 







[3.0 


-3.3] x 10~ 5 


D° 


-> pW 


9.7 x 10" 


-10 


[3.5 


-4.7] x 10~ 7 


D° 


— > LOfi+fT 


9.1 x 10" 


-10 


[3.3 


-4.5] x 10~ 7 


D° 









[6.5 


- 9.0] x 10- 8 


D + 


-» p W 


4.8 x 10~ 


-9 


[1.5 


- 1.8] x 10~ 6 


D+ 


K* W 


1.6 x 10~ 


-9 


[5.0 


-7.0] x 10~ 7 


D + 


K* + 







[3.1 


-3.7] x 10~ 8 


D° 


-> K*VV 







[4.4 


-5.1] x 10- 9 



Table 3. The D -> P£ + f~, (( = e, fi) branching ratios. 



d prr 



BR 



SD 



BR 



LD 



D u 


-> K u rr 







4.3 x 10- 


-7 




-> n + i + r 







6.1 x 10" 


-6 


D° 


-> n°e + r 


1.9 x 10" 


-9 


2.1 x 10- 


-7 


D° 


— > rjrr 


2.5 x 10- 


-10 


4.9 x 10 


-8 


D° 


j]'i + r 


9.7 x 10- 


-12 


2.4 x 10 


-8 


D + 


7T + tr 


9.4 x 10- 


-9 


1.0 x 10 


-6 


D s + 


K + rr 


9.0 x 10 


-10 


4.3 x 10 


-8 


D + 


K + rr 







7.1 x 10 


-9 


D° 


K°^ + r 







1.1 x 10" 


-9 



The experimental upper bounds for D — > Pl+l decay rates 
were recently improved by FOCUS with new upper 



bounds of 10~ 5 or less, close to the theoretical predictions [ 
I71 I16I . The allowed kinematic region for the dilepton mass 
ma in the D — > Pt + (~ decay is ma = [2mc, mo - flip]- The 
long distance contribution has resonant shape with poles at 
ma = m p a, m u , fl7l ll6l . There is no pole at ma = 
since the decay D — > Py is forbidden. The short distance 
contribution is rather flat. The spectra of D — > Pe + e~ and 
D — > P/i + fi~ decays in terms of ma are practically identi- 
cal. The difference in their rates due to the kinematic re- 
gion ma - [2m e ,2m M ] is small and we do not consider 
them separately. 

The differential distribution for D + >° n +Q t + r [ □ QH 
indicates that the high mass dilepton region might yield 
an opportunity for detecting c — > u€ + €~. The kinematics 
of the processes D — > W + £~ would be more favorable to 
probe the possible supersymmetric enhancement at small 
mu, but the long distance contributions in these channels 
are even more prominent [ 6 1 . 

We undertook an investigation of the D° — > yy decay 
and recently CLEO gave the first upper limit BR(D° — > 
yy) < 2.9 x 10" 5 [ 1171 . The short distance contribu- 
tion is expected to be rather small and therefore the main 
contribution is coming from long distance interactions. 
The total amplitude is dominated by terms proportional 
to a i that contribute only through loops with Goldstone 
bosons. At this order the amplitude receives an annihilation 
type contribution proportional to the ao Wilson coefficient, 
given by the Wess-Zumino anomalous term coupling light 
pseudoscalars to two photons. Terms which contain the 
anomalous electromagnetic coupling of the heavy quark 
Lagrangian are suppressed compared to the leading loop 
effects [0. The invariant amplitude for D° — > yy de- 
cay can be decomposed using gauge and Lorentz invari- 
ance into a parity conserving and parity violating parts [ 
|5J. Within this framework [ 8 J, the leading contributions 
are found to arise from the charged n and K mesons run- 
ning in the chiral loops. As a result of these internal ex- 
changes our calculation predicts that the D — > 2y decay is 
mostly a parity violating transition [[5). We find that the 
predicted branching ratio is 



BR(D -> yy) = (1.0 ± 0.5) x 10" 



(4) 



The short distance contribution within the SM has been es- 
timated in [ 16 1, which found that it cannot exceed the long 
distance part. 

Recently we have studied the dominant contributions to the 
rare decay mode D — > €*t~y [ 9 1. The long distance contri- 
butions can be divided into nonresonant and resonant con- 
tributions. The first is the result of n, K loops, while the 
second one occurs on account of the D° — > Vy — > £ + {~y 
transition where V = p°, w, n vector mesons. We have 
calculated that nonresonant contributions give the rates 



BR(D° 
BR(D° 



> ye + e-) nonres = 1.29 xl0-'°, 

> wV)nonres = 0.21 X 10" 10 . 



(5) 
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The resonant contributions are 



BR(D° -> py -> try) ^ 5 x 10~ 
BR(D° -> ajy -> try) ^ 8 x 10~ 
BR(D° -> cf>y -> fry) ^1x10 



-9 



(6) 



As seen from these results the largest contributions come 
from the <p resonance. Finally, combining both contribu- 
tions we estimate that the dominant long distance contri- 
butions in these decays give 



BR(D° -> ye + e~)sM = UX 10" 
BR(D°-> ry uV)sM = l.lXl0- 



10 



(7) 



Using simplest supersymmetric extension of SM we inves- 
tigate possible effects. The largest contribution is expected 
from gluino-squark exchanges rl9l l!2l . Using single mass 
insertion we find [|5J, 



BR(D° -» yte-) SM 
BR(D° 



1.4 x 10 



-9 



7A/V)sm= 1.2 X10" 9 . 



(8) 



D° — > K P with P = n , rj, rj' with nonresonant K n, which 
we consider to be the most likely candidates for early de- 
tection. Here we used the heavy quark chiral Lagrangian 
supplemented by light vector mesons, as the theoretical 
framework. 

The radiative K — > nny decays offer useful information in 
the study of the K meson decays. Their counterparts in 
charm sector are charm meson decays into P^y. Using 
the factorization approximation for the calculation of weak 
transition elements we use the information obtained from 
semileptonic decays N18II19I . The nonleptonic D — > Kjt 
amplitude cannot be calculated accurately in the factor- 
ization approximation from the diagrams provided by our 
model. Such a calculation gives a rather good result for the 
D + — > K°7T + channel but is less successful for the D° — > 
K + 7r~ decay. In order to overcome this deficiency and to be 
able to present accurately the bremsstrahlung component 
of the radiative transition, we shall use an alternative ap- 
proach. We take the experimental values for the D — > Kji 
amplitudes, in the calculation of the bremsstrahlung com- 
ponent. In order to accommodate this, we write [ 20 1 the 
decay amplitude as 



The situation rather changes if 7?-parity violation is al- 
lowed as in [ 9 1 . The effective interaction is then 



eff 



A 'i2k A 'ilk 

2M L 



(u L /c L )(4r^L), 



(9) 



where are the coefficients of the "R-parity breaking 
terms of the superpotential [ 9|. Using new experimen- 
tal bound BR(D + -> n + p?yT) < 8.8 x 10~ 6 [ 3 1 we obtain 
a new bound 



Ui < 0.003 



100 GeV 



(10) 



Using this bound and allowing photon emission from t 
and t, with photon energy cut E y > 50MeV and E y > 
100 MeV 



+ F 1 [(q-s)(p-k)-(p-s)(q-k)] 



(12) 



where Fq is the experimentally determined D — > K7r am- 
plitude and F\, F2 are the form factors of the electric and 
magnetic direct transitions which we calculate with our 
model. When intermediate states appear to be on the mass 
shell, we use Breit-Wigner formula. Thus, we get [ 1201 
for the branching ratios of the electric transitions including 
bremsstrahlung, with |Fo| determined experimentally and 
taking the photon energy cut E c > 100 MeV 



BR(D + -> RVy) 



PV.ex 



(2.3 - 2.5) x 10" 



(13) 



For the D° radiative decay we get 



-9 



BR(D° -» re + e-)? r£50MeV = 4.5 x 10 

BR(D° -» W V)fi ysnvuv = 50 X 10" 9 



v >50MeV 
y >50MeV 

BR(D° -» re + e-)^ £l00MeV 



-9 



BR(D° - 



= 4.5 x 10 
: 46 x 10~ 9 . 



(ID 



Note that the SM predictions are not affected by the cuts on 
the soft photon energy at the order E 7 > 100 MeV as the 
bulk of the contribution either comes from the resonances 
or the low p 2 region [ 9|. The cut on E y is the cut on the 
high p 2 region. 

In addition to the D meson rare decays in which the FCNC 
transitions might occur, we have studied the Cabibbo al- 
lowed radiative decays D + — > K°7r + y, D° — » K ~n + y and 



BR(D -> K -jr + y) 



PV.ex 



(4.3 - 6.0) x 10" 



(14) 



The uncertainty in the Fq/F\ phase is less of a prob- 
lem in D + -> K°7T + y than in D° -> KTn + y. If we 
take the bremsstrahlung amplitude alone as determined 
from the knowledge of |Fo|, disregarding the direct elec- 
tric F\ term, the above numbers are replaced by 2.3 x 10~ 4 
for D + decay and 5.5 x 10~ 4 for the D° decay. The 
contribution of the direct parity violating term (putting 
F = 0), is BR(D + -> K°7r + 7) dh ,p V = 1.0 x 10~ 5 and 
BR(D° -> K ^ + y) dir ,p V = 1.64 x 1(T 4 . For the parity con- 
serving direct magnetic transition we get [ 20 1 

BR(D + -> KVy) PC = 2.0 x 10~ 5 , (15) 
BR(D° -> K"7r + y)pc = 1.4 X 10~ 4 . (16) 
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Hence, the two direct transitions are predicted to be of 
about the same strength. If we disregard the contribution 
of vector mesons to the direct part of the radiative decays, 
the parity-conserving part of the amplitude is consider- 
ably decreased, by one order of magnitude in the rate in 
D + — > K°7t + y decay and by two orders of magnitude in 
D° — > K ~n + y. On the other hand, their contribution is not 
felt in a significant way in the parity - violating part of the 
amplitudes. In any case, the detection of the direct part of 
these decays at the predicted rates, will constitute a proof 
of the important role of the light vector mesons. 

The decay mode D° with neutral pseudoscalar mesons has 
no bremsstrahlung part of the amplitude. The light virtual 
vector mesons give the main contribution to the decay am- 
plitude. The decay D° — > K°ny is predicted [ 20 1 to have 
the rate 

BR(D° -> K°7ry) = 3.4 x 10~ 4 , (17) 

with comparable contributions from the parity-conserving 
and parity-violating parts. The decays D° — > K°j](ri')y are 
expected with rates 1 . 1 x 10 and 0.4 x 10~ 7 , respectively, 
and are mainly parity conserving. 

In summary: we find that the best candidate to search for 
new physics effects in charm decay seems to be the dif- 
ference BR(D° p°y) - BR(D° -> uy). The observation 
of the rate of BR(D° — > yy) larger than 10~ 8 can signal 
new physics, while the rate BR(D° — > fi + fi~y) of the order 
10 -8 can be explained as a result of the 7?-parity violating 
MSSM effects. We hope that interesting features of these 
decays result in an experimental search in the near future. 
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